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We investigate a model for the regularly depleted two-leg 
spin ladder systems. By using Lieb-Schultz-Mattis theorem, 
it is rigorously shown that this model realizes massless exci- 
tations or, alternatively, a degenerate ground state, although 
the original spin ladder system has a spin gap and a unique 
ground state. The ground state of the depleted model is ei- 
ther a spin singlet or partially ferromagnetic reflecting topo- 
logical properties of the depleted sites. In order to show that 
the low-energy excitations are indeed massless, we proceed 
our analysis in two different ways by resorting to effective 
field theories. We first investigate an effective weak-coupling 
model in terms of renormalization group methods. Although 
the tendency to massless spin excitations is seen in the strong 
coupling regime, it turns out that the model is still massive 
for any finite coupling, implying that a conventional weak- 
coupling approach is not efficient to describe massless modes 
in our model. To overcome this difficulty, we further study 
low-energy properties of the depleted spin model by mapping 
on the non-linear sigma model, and confirm that the massless 
spin excitation indeed occur. 

75.10.Jm,75.10.-b 



I. INTRODUCTION 

Since many decades one-dimensional quantum spin 
models have been a prominent subject of theoretical 
study, representing a simple, but very rich class of many- 
body systems. During the last three years the problem of 
several coupled quantum spin chains has attracted much 
attentionE One reason lies certainly in the discovery of 
quasi one-dimensional spin systems in layered cuprate 
materials which are the basis of high-temperature super- 
conductors. Synthesized under high pressure the Cu02- 
planes of the infinite layer compound SrCu02 creates 
regular line defects which lead to nearly decoupled one- 
dimensional spin systems with ladder structure. There is 
a homologous series of compounds, Sr n _i Cu„+i O211, 
which form ladders with -different numbers of legs (= 
(n+l)/2) depending on n.H The only (two-leg) spirulad- 
der compound previously known was (VO)2 P2 O7.B 

Spin 1/2 ladder systems with antiferromagnetic (AF) 
nearest neighbor coupling behave differently if they have 



an even or an odd number of ladder legs. In case of 
an even number of legs the system has a resonating va- 
lence bond (RVB) ground state (short-range singlet cor- 
relation) and a gap to the lowest excitations. On the 
other hand, an odd number of legs leads to gapless ex- 
citations and a ground state with quasi long-range order 
similar to that of the single AF spin 1/2 chainjj These 
properties have been clearly observed in experiments for 
the compounds SrCu203 (two-leg ladder) and Sr 2 Cu305 
(three-leg ladder) and the size of the excitation gap for 
the two-leg ladder system is in good agreement with the- 
oretical predictions.B 

An interesting new aspect occurs if these systems con- 
tain non-magnetic impurities, i.e., some spins are re- 
moved from the ladder. This can be achieved, for ex- 
ample, by substituting Cu-ions by non-magnetic Zn. Ex- 
periments showed that even a rather small concentration 
of Zn (> 1%) is sufficient to yield a transition to antifer- 
romagnetic long-range order in this two-leg ladder com- 
pound SrCu203 as observed in recent experiments!] In 
model calculations it was actually demonstrated that in 
the vicinity of a removed spin a staggered magnetization, 
develops which extends over several lattice constants.Q 
Each removed spin leaves locally an unpaired spin 1/2 
degree of freedom behind. These residual spins interact 
with each other through exchange of excitations of the 
ladder system analogous to the RKKY-interaction. It 
can be shown easily that if two impurities lie on the same 
sublattice (ladders are bipartite lattices) then the effec- 
tive interaction between the corresponding residual spins 
is ferromagnetic (FM). On the other hand the interactio 
is AF if the two impurities occupy different sublattices. 
In the ground state these residual impurity spins corre- 
late due to their effective interaction. The orientation of 
the local staggered moments associated with each impu- 
rity is also correlated with that of the impurity spin. One 
can easily see that this leads to an in-phase alignment of 
the local staggered moments throughout the whole sam- 
ple. In this context it is crucial that this system has no 
frustrating interactions. With coupling among the lad- 
ders this behavior is indeed sufficient to create AF long- 
range order. For the coherence of the AF correlation it 
does not matter whether the impurities are located in a 
regular or random way along the ladder. However, it is 
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still not clear yet whether and under which conditions 
true long-range order wuld also emerge in the ground 
state of a single ladder BeI 

While the pure two-leg spin ladder has an excitation 
gap, a finite impurity concentration seems to yield gap- 
less excitations. This was argued recently for the random 
depletion of spins, where for low energies the system can 
be reduced to an effective random spin chain.Ej More re- 
cently, lino et al.liiJ and Motome et ahcll analyzed the 
excitation spectrum numerically. For very small impu- 
rity concentrations the spin gap feature of the two-leg 
ladder system is still a dominant structure and only few 
low-lying excitations appear far below the gap. However, 
already at a concentration as small as 4% a crossover to 
a regime occurs where-the original excitation gap has es- 
sentially disappeared.EJO 

In this paper we would like to consider the problem of 
low-lying excitations on a more rigorous basis for the case 
of the regular arrangement of impurities in two-leg ladder 
systems. In sec II, we first prove rigorously by the Lieb- 
Schultz-Mattis theorem that any finite concentration of 
removed spins leads to an excited state with the energy 
of 0(1/ N), which implies the formation of either gap- 
less excitations or the existence of a degenerate ground 
state with spin gap (N: number of lattice sites of the 
finite system). In the next step, we will attempt to de- 
termine which possibility is actually realized in the model 
by field-theoretical methods. For this purpose we first in- 
troduce in sec. Ill an effective Hubbard-type ladder model 
and study its low-energy properties using the renormal- 
ization group method. In this Hubbard model the ef- 
fective depletion of spins is incorporated in the strong 
coupling limit of on-site interaction which is introduced 
to produce the singlet bound states at each depleted lat- 
tice point. Though the tendency to gapless excitations is 
seen as the singlet bound state becomes strong, it turns 
out that the model is still massive for any finite couplings. 
This may suggest that it is not straightforward to incor- 
porate the effects of depletion in a conventional contin- 
uum limit. To clarify this point, we further study in sec. 
IV the depleted spin model by employing a complemen- 
tary field theoretical approach, i.e. the non-linear sigma 
model. To this end, we first calculate the dispersion of 
the spin waves by the Holstein-Primakoff mapping. We 
then find that in the case of singlet (partially ferromag- 
netic) ground state, the spectrum of the lowest band is 
linear (quadratic) without a gap, suggesting that correla- 
tions between unpaired spins are in fact essential for low- 
energy excitations. Based on this observation, we study 
the model with singlet ground state by using the map- 
ping to the non-linear sigma model. It is found that the 
coefficient of the topological term in the effective sigma 
model is ni, analogous to the spin 1/2 antiferromagnetic 
Heisenberg chain, suggesting that our spin ladder sys- 
tems with periodic depletion indeed have massless spin 
excitations, being consistent with the analyses in sec. II 
and III. Brief summary is given in sec. V. 

We wish to mention here that the model studied nu- 



merically by lino and Imadalld is similar to the one we in- 
troduce in Sec. II. They discussed how the massless states 
appear as the impurity concentration is increased. Our 
approach based on the rigorous statements may provide 
the results complementary to theirs, and also via the 
present analysis we will clearly see how the coherence 
of the eigenfunction is modified by the depletion. 

II. REGULARLY DEPLETED SPIN LADDER 
A. Model 

As mentioned above, the presence of the gapless stateis. 
a quite common feature in the depleted spin ladders.Qtil 
Although this has been already claimed by various stud- 
ies based on impurity models, it may be important to 
rigorously prove that such a gapless state can be indeed 
realized by the depletion. In order to study the effects of 
depleted spins on two-leg ladder systems, we study here 
a special class of the spin models, i.e., regularly depleted 
spin ladder systems, from which one can clearly see the 
drastic change of the ground state. The model Hamilto- 
nian we consider is 

H = + if (2) + H coup , 

N-l 

HV=J^{l-W){l-1$Js¥>-S$ lt 

3=0 
N-l 

3=0 

where 9 is defined by O^p = 6j and — 9j +m /2 with 

q . _ f 1 for j = (mod m) ,^ ^ 

3 I otherwise ' \ ■ ) 

In what follows, we set J\\ — 1 and Jj_ — J for simplicity. 
This model is defined on the two coupled chains, each of 
which consists of N sites labeled by j = 0, 1, • • • , N — 1. 
An even-integer m denotes the period of depleted sites 
(we call them impurities) with N = mM , where M de- 
notes the number of impurities on each chain. Periodic 
boundary conditions are imposed on the system. We 
show the model in Fig. |l| schematically for the case of 
m = 6. 

Although our model seems rather special, it is to be 
noted that the model contains an essential property ex- 
pected generally for the depleted models£2l: the gapless 
states are essentially formed by the coherent motion of 
unpaired spins generated by the depletion,Q though an 
unpaired spin may be a rather complicated object in gen- 
eral. For example, in the case where the period of the 
depletion is large (dilute limit of impurities), our model 
reproduces the known results for the system with a few 
non-magnetic impurities. But also for the high density 
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FIG. 1. Regularly-depleted two-leg spin ladder system with 
m = 6. The number assigned to each site is used in the 
spin-wave analysis in section IV. A. 



limit, we can discuss what is essential for the formation of 
a gapless state in the depleted systems. Moreover, since 
we can give a rigorous statement for the model, it may 
serve as a key model by which we can check whether an 
approximate treatment works well or not when we study 
wider class of depleted ladder systems. In fact, we will 
treat a slightly different model in the next section, for 
which the exact results obtained in this section provide 
a guideline for analyzing the results correctly. 



B. The Lieb, Schultz and Mattis Theorem 

By using the Lieb, Schultz and Mattis (LSM) 
theoremtd, we now prove that the depleted model (2.1) 
indeed has an excited state of 0(1 /N). The LSM the- 
orem is elementary, but provides quite powerful tool to 
demonstrate that a given system is gapless or has degen- 
erate ground states. This theorem i-Ga&i be applied not 
only to half-odd integer spin modelaiStid but also to the 
spin ladder systems, which suggests that odd-leg ladder 
models are expected to have massless excitations, while 
it is not applicable to even-leg ladder models, indicating 
the presences of an excitation gap. The relation between 
the LSM scheme and the massive RVB states was dis- 
cussed previously in detail in literature (see for example 
Ref.14 and references therein). 

We n ow wish to apply this theorem to the depleted 
model 



(2.1 



The theorem uses two properties. One is 
the symmetry of the model based on the translation com- 
bined with reflection. The other is the response of the 
ground state energy to twisted boundary conditions, or 
equivalently to external gauge fields. This plays an im- 
portant role for studying low-energy properties in quan- 
tum many-body systems. In particular, such effects ap- 
pear in quite different ways for massless and massive sys- 
tems. For example, for massless systems, the ground 
state energy increases as a function of the twist angle of 
0(1/ N), where N is the size of the system, so that the 
ground state changes into one of the excited states when 
we follow the spectaji-43ow of the ground state up to the 
twist angle $ = 27rt3'll3. On the other hand, for massive 
systems, there is a gap above the ground state of 0(1), 
and hence we cannot reach the excited states by twist- 
ing the ground state. In this way, we can use twisted 
boundary conditions to find a gapless excitation. 



We now turn to the question whether we can find such 
an excitation of 0(1/ N) in the present model. To this 
end, let us introduce the operator T TO / 2 for m/2-latticc 
translation and the reflection operator R, 



T 



rn / 2 j m/2 — KJ j+m/2> 

RS^Rr l = S^. 



(2.3) 



Let O be the product operator O — RT m / 2 - The total 
Hamiltonian is invariant under the operation of O since 
0# (i,2) -i = #(2,1) and OH^pO- 1 = ff coup . On the 
other hand, for a finite size system the ground state is 
non degenerate within the subspace of fixed S* ot (= 0), 
so that we can write O|$o) = e* a |$o)- 

Now we introduce the twist operator, which plays the 
central role in the LSM theorem, 



exp 



2m 



N-l 
3=0 



(2.4) 



This operator transforms under O such that 



OU^O- 1 = U {2) exp 



l£(2)* 



x exp 



N 

m/2-l 

2m £ 



(2) W (2) 2 
3 



0[ /(2) -i = [/(i) exp 



3=0 

raui 



x exp 



N 

m/2-l 

2m E (i~0f } )S 



j 



3=0 



(2.5) 



and therefore U — U^'U^ is transformed as 



, / mm 
QUO- 1 = Ucxp ( —S' 



x exp 



-U. 



N 

2 m/2-1 
i=l j=0 



3 



(2.6) 



The second line comes from the facts that we are con- 
cerned with the total spin S z = sector and that in the 
last exponential there are m — 1 (odd) spins inside of the 
dashed box in Fig.[iJ each of which gives the factor — 1. 
Owing to this property, we can construct an excited state 
| VP) = U\^o) with different "parity", which is orthogonal 
to the ground state, (\P |*) = 0. 

The remaining task is to calculate the energy incre- 
ment SE = - (*o|-ff|#o}< We have 

2 N-l 

5£ = ^^(l-f)(l-ft) 

z=l 3=0 



3 



x(cos^-l) (*„|S<" + Sf-|*.) 

<2(A--M,(l-c„ s f)< 2 (l-|)£. (2.7) 

Namely, the excited state orthogonal to the ground state 
has the excitation energy of 0(1/N). Therefore, we end 
up with the rigorous statement that our model ( |2 . 1[ ) for 
regularly depleted ladder systems has either the massless 
spin excitations or the ground state is degenerate in the 
thermodynamic limit. In subsequent sections, we will 
focus our attention to the question whether low-energy 
excitations of our system are indeed massless, and, if so, 
which type of spin excitations contributes to this massless 
mode. 

Before concluding this section, we wish to remark the 
following points. As mentioned above, it is the most cru- 
cial point in the LSM theorem whether we can construct 
a state orthogonal to the ground state by the use of the 
twist operator. In the present system, an unpaired spin 
in the dashed box in Fig.[j] gives the factor —1. In the 
ordinary two-leg ladder systems without depletion, we 
have always an even number of spins inside the dashed 
box (however, in this case translation by one site is suffi- 
cient to consider, which still gives the factor (— l) 2 = 1). 
Therefore, we can explicitly see that a drastic change in 
the ground state directly reflects how the phase coherence 
of the wave function is changed by the depletion. 

Some aspects of the ground state can be discussed 
by applying Marshall's theorem to our modelOiia The 
ground state has the spin quantum number 5 = (M) 
for m = 4 (2) mod 4, so that the ground state of our 
model is either a spin singlet or partially ferromagnetic. 
Therefore, in the partially ferromagnetic case, the ground 
state itself is degenerate, since the state under consider- 
ation (S z = state) is, of course, a member of the spin 
S = M multiplet. As mentioned in the introduction, this 
can be easily understood by noticing that the effective 
interaction between the unpaired spins is ferromagnetic 
(antiferromagnetic) if the depleted sites are on the same 
(different) sublattice(s).El 

III. WEAK-COUPLING APPROACH TO 
HUBBARD LADDER 

In the previous section, we have proved that our model 
for the regularly depleted ladder systems is characterized 
by massless spin excitations or degenerate ground states 
with a gap. In order to determine which possibility is ac- 
tually realized and to study, if possible, the nature of ele- 
mentary excitations in more detail, we wish to introduce 
a low-energy effective theory in the continuum limit and 
study its low-energy properties by using the renormal- 
ization group (RG) method. One immediately notices, 
however, that it is not straightforward to study the effect 
of the periodic depletion in an ordinary bosonization ap- 
proach. For example, we cannot take a conventional way 



in which each chain is first bosonized and the inter-chain 
coupling is then taken into account via the RG procedure, 
because a bare single chain would be divided into many 
disconnected pieces by the depleted sites. To avoid this 
difficulty, it is desirable to find a way to include the effect 
of the periodic depletion after bosonization. Based on 
these observations we introduce a Hubbard-type ladder 
model in which the effect of the depletion is incorporated 
in terms of the on-site interaction. We then investigate 
its low-energy properties by the one-loop RG method. 
We wilkuse the notations similar to those of Balants and 
Fishciila, who studied the ordinary Hubbard type ladder 
systems without depletion. 

A. Model 

The Hamiltonian we consider is the two-leg ladder of 
correlated electrons, 

2 

H = ^2( H hop + H int) + #coup, 

i=l 

3=0 T 

JV-1 
3=0 

jV-1 

tfcoup = J$>- - 0f ] )sf ] ■ Sf\ (3.1) 

3=0 

where we have introduced two types of the Hubbard in- 
teractions U > and U + U' < 0. The interaction U is in- 
troduced to produce Heiscnberg spins on the lattice sites 
for large U . On the other hand, the attractive interac- 
tion U' is introduced to effectively represent the depleted 
spins. Namely, in the U + U' — > — oo limit, electrons on 
j = (rn/2) mod m sites for chain-1 (-2) form bound 
states, i.e. onsite spin singlets, and hence the spin de- 
grees of freedom are effectively frozen out on these sites, 
as illustrated in Fig.||. In order to correctly reproduce the 
depleted ladder systems, the number of electrons should 
be 2(N + M), and then the corresponding Fermi momen- 
tum is 

h * = i( 1 + ^)' < 3 - 2 > 

where clq is the lattice constant. 

We now investigate low-energy properties of this model 
with the RG method. Passing to the continuum limit, the 
fermion operators become 

c$/V^ - e- lkFX 4>iL*(x) + e lkFX ^ lRa {x) (3.3) 
with aj — > x, and the hopping Hamiltonian is reduced to 
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FIG. 2. Schematic illustration of the two-leg Hubbard lad- 
der in the strong-coupling limit U' — > — oo 



{JihJiL + JiRJiR) 



7 



(3.4) 



where up = 2tao cos(7t/2to). Here and in what follows, 
operator products should be normal-ordered, though we 
do not explicitly indicate this. J^ L r R s is the left (right) 
component of the U(l) and SU(2) currents defined by 



(3.5) 



where tx° is the unit matrix: The component /x = is 
the usual U(l) current divided by 2, and those with a = 
fi = 1,2,3 are the SU(2) currents, respectively. Next we 
bosonize the Hubbard interaction in this approximation. 
In the following we neglect the oscillating terms which 

(i) 

are incommensurate with the Fermi momentum and #) , 
as they may disappear after the integration. According 
to calculations briefly summarized in Appendix A, the 
intra-chain interaction leads to the following Hamiltonian 
density, 



where 



"ftint — 9pJihJiR + 9x P {J\lJ2R + J^L^Ir) 

+g^(-r\M tL M} R + M} L M lR ) 

+9crJiLJiRi 



(3.6) 



(3.7) 



and we have dropped the chiral interactions because they 
merely renormalize the Fermi velocities. Here, repeated 
indices i and a are summed over. Note that the Umklapp 
interaction appears though the density of the system un- 
der consideration is more than half-filling. This is due to 
the existence of 9j in eq.(3.1), which fixes U' on regularly 
placed lattice sites. We can also take the continuum limit 
of the inter-chain coupling 

— Hcoup = 9xa{J\hJ2R, + J^hJlh) + dt^ijAJ^R + L L ^r) 

+g tp {LlL^+L^L R ) (3.8) 

where ii(_R) is the staggered components of the spin op- 
erator 



re 

n L{R) 



VlL(R) a ^rV2L{R)f3 



(3.9) 



The initial values of these coupling constants are given by 
the Hubbard and inter-chain interactions, as summarized 
in Table J, where we have defined U' = U' /m, U = U + 
U' and J = J(l — 2/m). The bosonized forms of the 



9p 9%p 


9u 


g<? 


9xa 


9ta 




-2Ua 


—U'ao 


2Ua 


—Jao 


— 2 J 0-0 


§ Jao 



TABLE I. Coupling constants 
interactions are also listed in Appendix C. 

B. RG Equations and Flows 

We now derive the RG equations for the above cou- 
pling constants in one-loop order. Using the operator 
product expansions and resultant RG equations given in 
Appendix B, we end up with the following set of scaling 
equations, 



dg P 

dl 

dg xp 

dl 
dg u 

dl 
dg G 

dl 

dg X(y 

dl 
dg t g 

dl 



dgtp 
dl 



-4g 2 u-\(9l + 3gl), 
■-;(9tp + 3 9ta), 



= -g P 9u, 



;{9tc + 9tp)9tc 



= ~9xa - 2^ta ~ 9tp)gt. 



:(g P - gx P ) + (g* + gxa) 



9to 



-^{g<y- 9xa)gt P , 

~2^P ~ 9xp)gtp - ^(.9<x - 9xa)gta 



(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 

(3.15) 
(3.16) 



with I — InL, where gu = 'gk/^VF- 

We have numerically integrated the above set of equa- 
tions to obtain the RG flows. Typical examples of the 
RG flows are shown in Figs.|^ and [l| We note that a 
characteristic scale l c exists, at which all renormalized 
couplings except for g xp exhibit divergence properties, 
though those in Fig.[| look convergent at a first glance. 
We will see below that this divergence is mainly driven 
by the charge sector, so that l c may be identified with the 
scale at which the mass gap is generated for the charge 
sector. Accordingly, couplings for other spin modes are 
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FIG. 3. RG flows for (a) g p , (b) g u , (c) g t<7 and (d) g tp : 
Plotted values are multiplied by the factor 10. g xp ~ in this 
figure. The parameters used are U = 1, U' = — 20, j = 0.01 
and m = 50. 




FIG. 4. RG flows for (a) g a and (b) g xa : Plotted values 
are multiplied by the factor 400. The parameters used are 
the same as those in Fig|j. 

strongly affected. In what follows, we analyze character- 
istic properties of RG flows around l c analytically. Note 
first that the key parameter in our approach is U', which 
would make singlet bound pairs on particular sites, freez- 
ing out the spin degrees of freedom. Now we fix the other 
parameters, e.g. U — 1, m = 50, J = 0.01. In this case, 
the Umklapp interaction turns out to be relevant in the 
range U' c < V < ( U' c ~ —40 for this particular choice 
of interactions. We will explain the meaning of U' c be- 
low). In this ran ge, g is most relevant and we can neglect 
g t<J and g tp eqs.( 3.10 )~( 3.12|) , resulting in 



dg P 

dl 
dg u 
dl 



-9 P 9u, 



(3.17) 



near the divergent point. This is a familiar U(l) scaling 
equation, from which we obtain g p ~ —l/(l c — I). This 
expression near the divergent point is confirmed by the 
numerical integration of the RG equations. This result 
means that the charge gap is generated by the Coulomb 
interactions U and U' . 

What happens for the spin sector in this case ? Note 
that the next relevant interactions are the inter-chain 



staggered interactions g tC r and g tp for the spin sector, 
since they act as composite of spin and charge degrees 
of freedom, the dimension of which is unity if we replace 
the massive charge degree of freedom by its mean value, 
as can be seen in Table|n]. Neglecting small terms, we 
thus obtain 



dgta 
dl 

dg tp 
dl 



1 

1 

-^9 P 9t P - 



(3.18) 



Using the above-obtained g p for the charge sector, we 
have the RG flows 



9ta 

9t P 



9ta0 



(lc-1) 1 / 2 ' 

(ic-iy/ 2 ' 



(3.19) 



near the divergent point. Therefore, it is seen that the 
spin modes scale to the strong coupling fixed-point and 
the spin gap is still open. This may not be plausible for 
the depleted spin ladders, according to the rigorous theo- 
rem obtained in the previous section. This is because the 
effect of the depletion is not fully.-taken into account in 
the present one-loop RG method.Ej It should be noted, 
however, that the effect of the depletion actually has the 
tendency to make the spin sector massless. Namely, we 
can indeed see that the coefficients gt&o and gt P Q quickly 
decrease with increasing |J7'|£3 In fact from the numer- 
ical integration of the RG equations, we estimate them 
as 



9ta0, gt P o 



const. 
\U'\ a 



(3.20) 



with a ~ 1 (See Fig.g). Therefore, although the spin gap 
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FIG. 5. (a) —gto-o X 10 3 and (b) g tp o x 2 ■ 10 2 as functions 
of negative U' . 

still exists for finite U', we can expect, by smooth extrap- 
olation to U' — > — oo, that the divergence in the above 
RG flows may be suppressed, leading to the massless spin 
modes. 

We should mention here that the above treatment 
would fail at some negative critical value V = U' c if 



G 



we only increase \U'\ by keeping U fixed. Namely, for 
\U'\ > U' c , the charge mode becomes massless, and then 
we cannot use the model as a weak-coupling analog of the 
depleted lattice. This pathological result may come from 
the fact that the extremely strong attractive U' may over- 
whelm the repulsive U completely, resulting in the spin 
bound states similarly to the attractive Hubbard model 
and hence making the charge mode massless. So, it may 
be necessary to take a strong-coupling limit by appropri- 
ately tuning the ratio of U and U' in the present model 
with attractive U' . 



C. Discussions 

Within the above analysis for the effective weak cou- 
pling model, it turns out that the model is still massive 
for any finite coupling, though the tendency to gapless 
excitations is seen as the singlet bound state becomes 
strong. This may suggest that it is not straightforward 
to incorporate the effects of depletion in a conventional 
continuum limit. In this sense, we should say that our 
analysis in this section is not satisfactory to describe 
massless excitations in the depleted spin model. In this 
connection, we would like to note here that the present 
RG analysis is based on the effective weak-coupling "elec- 
tron" model because we believe that this electron model 
belongs to the same universality class of the original spin 
model. In fact we have numerically solved the RG equa- 
tions not only for the spin sector but also for the charge 
sector, and found that all the trajectories flow to the 
strong coupling fixed points. Furthermore, as an exten- 
sion, one could examine whether the above situation is 
improved if we solve the RG equations only for spin sec- 
tors by discarding the charge degrees of freedom, which 
should be a more accurate approximation for the original 
spin model. Unfortunately, even if we concentrate on the 
RG equations for the spin sector, we still encounter the 
problem that the spin sector is massive since the dimen- 
sion of the inter-chain couplings reduce to 1, as can be 
seen from Table [n]. In this sense, a conventional weak- 
coupling bosonization approach, in which two chains are 
assumed to be weakly coupled with each other, does 
probably not allow a straightforward description of the 
low-energy properties of the depleted model. We thus 
have to find or develop a more effective way to incorpo- 
rate the effect of depletion correctly. This problem will 
be discussed in the following section. 

IV. NON-LINEAR SIGMA MODEL APPROACH 

In the previous section, we have seen that it is not 
straightforward to incorporate the effect of depletion in 
an ordinary weak-coupling approach. In particular, it 
is quite difficult in this weak-coupling model to figure 
out which kind of spin excitations would actually become 



massless. To clarify this point, and also to confirm that 
the massless modes are actually realized, we would like 
to consider a different field theoretical .anproacb. by using 
non-linear sigma model techniques B3~Eji9 

Let us reconsider here the essential ingredients impor- 
tant for low energy excitations in our model. When ap- 
plying the LSM theorem in sec. II, we have seen that 
the presence of unpaired spins associated with vacant 
sites play a crucial role. All other spins are essentially 
bound into local singlet pairs and cannot contribute to 
the very low energy spectrum. The vacant sites break the 
uniformity of the system by introducing these unpaired 
spin degrees of freedom and by the local polarization of 
a staggered moment in their vicinity. A short-ranged 
effective interaction among the unpaired spins is medi- 
ated through polarization of the remaining spins, which 
is (anti-)ferromagnetic for m = 2 (4) mod 4. These inter- 
actions are weak and are assumed to introduce a gapless 
spectrum. In a weak coupling approach in the previous 
section, we have started with ordinary boson fields by 
simply applying the conventional bosonization schemes 
to each chain, and have not explicitly taken into account 
the above properties in the beginning. Therefore, in order 
to describe the massless mode better, we need a rather 
complicated combination of these boson fields by incor- 
porating higher-order interactions. This may be a reason 
why it was difficult in the weak coupling model to explic- 
itly construct the massless spin excitations. 

Based on these observations, we would like to consider 
below a complementary approach which allows us to de- 
scribe the formation of massless excitations more easily. 
We will concentrate on the regime of comparatively small 
to where the background of the staggered magnetization 
plays an important role. We first perform a spin-wave 
analysis of the depleted model to clarify the nature of 
spin excitations. We then convert the lowest spin mode 
to the non-linear sigma model, finding that the massless 
excitations are indeed generated for the depleted model. 

A. Spin- Wave Analysis 

We consider now a regularly depleted spin ladder with 
rather small m. Then the polarization cloud of staggered 
magnetization around each vacancy have a strong over- 
lap with each other. Assuming coherent, though slightly 
inhomogeneous staggered magnetization throughout the 
whole ladder we may use the Holstein-Primakoff map- 
ping. Because after regular depletion (period to) the 
unit cell of the ladder is large we represent the spin 

(i) 

operators by 2m kind of bosons a ; with i = 1 , 2 and 
I = 0,1, ■■•to — 1 (We here include the depleted sites 
i = 1, 1 = and i = 2, I = to/2 for simplicity. But their 
sector easily decouples, which has always zero-energy), 

^ +i =±Q-n«(mi + 0), (4-1) 
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1 — n 



■l)af\mj + l) 



(mj 

(mj +l)Jl- nf* (mj + I) 



(4.2) 



where ni is the number operator. In (4.2) we take the 
upper (lower) relation for i + I = odd (even). In the 
momentum representation, 



a, \mj + I) 



1 



M-l 



ex p 



fc=0 



2ni 



(mj + l)k 



a[ l) (k), 



(4.3) 



the Hamiltonian becomes, up to quadratic order of the 
boson operators, 



2m 
k i,j=l 

(4.4) 

where we have omitted some trivial constant terms. We 
have used a simplified notation (See Fig.jl] for this num- 



bering), d£\k) = ai, ■ 

(2) flA 

0-m+l, ' "j ' ' ') Cl m _ 1 {K) 



,(!) 



l(-fc) 



,( 2 ) 



a2m, and the 2m x 2m matri- 
ces h and A are given in Appendix D. By the Bogoliubov 
transformation 



2m 

E 

3=1 



which should satisfy the relations 

uv f -vu = o, uu f - vv 1 



1, 



(4.5) 



(4.6) 



we can reach the diagonal form of the Hamiltonian, up 
to constant terms, 



H 



3=1 



'jbjbj + const. 



(4.7) 



In Figs.|^J^, we show the dispersion uji as a function of 
k, in which we omitted the zero-energy modes associated 
with depleted spins. We can observe that several bands 
appear because of the periodic depletion. It is remark- 
able that within this approach the difference between the 
cases m — 2 and m = 4 (mod 4) is reproduced clearly. 
As discussed in section II, the ground state of the former 
case is ferromagnetic (partially spin polarized), which is 
reflected in the quadratic dispersion of the lowest band 
for the small momenta (Fig.||). In contrast the latter 
behaves quite differently, displaying a linear low-energy 
spectrum. Therefore, the lowest band is indeed due to 
the interaction between unpaired spins, as discussed at 
the beginning of this section. Though the lowest spin 
mode does not have an excitation gap in both cases, one 
cannot naively conclude a gapless spectrum for the sys- 
tem, because higher-order quantum corrections play an 




100 



FIG. 6. Dispersion of m = 6 system with J = 1 as functions 
of k with N — 600. In the usual notation, k — 100 point 
cor resp onds to the momentum 2ivi x 100/600 = ir/3. See 
eq.dDq) 
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FIG. 7. Dispersion of m = 8 system with j = 1 as functions 
of k with N = 800. In the usual notation, k — 100 point 
corresponds to the momentum 2ni x 100/800 = 7r/4. 



important role to determine whether or not the system 
is indeed gapless. This problem is investigated in the 
subsection B. 

We would like to mention here that although our ap- 
proach is valid for the case of small m, the qualitative 
properties of the low-energy spin excitation spectrum are 
the same for large m. In this case, however, a quantita- 
tive description would have to include the effects of the 
coupling among elementary excitations beyond (4.4) and 
(4.7). It is not our aim to address this problem here. 



B. Mapping to the Non-Linear Sigma Model 

Based on the above spin-wave analysis, we would like 
to now consider the lowest band due to the unpaired spins 
and to discuss the role of quantum effects by employing 
sigma model techniques. The results in the previous sub- 
section suggest that the system with m = 4 (mod 4) can 
be mapped to the non-linear sigma model, as is usually 
the case for uniform spin chains as well as uniform spin- 
ladder systems. We restrict ourselves to such cases in 
this section. 
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We use here the coherent state path-integral method.c3 
The partition function Z = tr exp(— 0H) can be rep- 
resented by Z = JH^ViniN^ij^expi-S), where 

the unit vector N^'(J) specifies the coherent state of 
the spin at the jth site of the ith chain such that 
(N^(j)\sf ] \N^(j)) = sN®(j) with s = 1/2, and 
fi(N) is the invariant measure J d(i(N)\N)(N\ = 1. The 
Euclidean action S is 



S = Sb + S\\ + S±, 



(4.8) 



where 



2 JV-1 



»=1 j=0 
r /3 2 



5n = ,s 2 



t=l J=0 



JV-1 



xiV«(i)-JVW(i + l) 



5j_ = Js 2 I drJ2(i-ef)(i-e^ 

Here w[7V] is the Berry phase defined by 

w[N] = fdr f duN ■ (d T N x d u N). 
J o Jo 

After staggering the configuration 

^%') = (-) i+i+1 ™ w (i), 

we suppose 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



In what follows, we assume the fields m and I are suffi- 
ciently smooth functions. If we want to describe not only 
the lowest band but also the massive bands, we must in- 
clude several kinds of fields. However, we concentrate 
on the lowest one only. We expect that the above (non- 
trivial) assumption are sufficient for this description. 

Let us now derive the continuum limit of the action. 
First we consider the Berry phase term and divide it into 
two contributions, 



2 JV-1 

-* s EEh i+j+ m™ w (j)] 

i=l j=0 
2 JV-1 

z=l j=0 



(4.13) 
(4.14) 



Taking now the continuum limit in each chain in 
cq.( 4.13| ), the topological terms cancel out due to the 
factor (— y, and we find 



(4.13) = 2is Jd 2 xl ■ (m x d 2 m), 



(4.15) 



where x\ = x (= aoj) and x 2 = t. Then the contribution 
from cq. ( 4.14 ) is evaluated 



M-l 



(4.14) =is [u[n (1) (mj)} - uj[n {2) (mj + m/2)] | 

M-l - 

~ is dT(—)Sn( 1 \mj) 



3 = 1 



n 



(1) (mj) x d T n {1) {mj)] , (4.16) 



where Sn^ (mj) 



n^(mj) — n^ 2 \mj 



— ^aQd x m{mj) + 2aol(mj) + O(oq). We have assumed 
that the short-range AF order is strong enough to justify 
our continuum limit. As noted above in this subsection, 
we here concentrate on the case m — 4 (mod 4), because 
the above equation is valid only for m/2 = even. Substi- 
tuting this relation, we have 

(4.14) = ~ (d 2 xm ■ (dim x d 2 m) 



2is 



d 2 xl ■ (m x d 2 m). 



(4.17) 



Next, we evaluate the intra-chain interaction. For ex- 
ample, the contribution from chain 1 is 



M-l m-2 



^ = Y r T E E [ nil) ( m J + - n^(mj + 1 + 1 
J 3=0 1=1 

2 m— 1 „j „ 

~^-J2J^ jdr[(d x m) 2 -A(-) l d x m-l + Al 2 ] 



QqS m 



d 2 x [(d x m) 2 + 4l 2 ] 



(4.18) 



This type of approximation may be justified for rather 
small m where we can assume that the region between 
the vacant sites is clearly dominated by the staggered 
spin polarization. Then all the intermediate spins vary 
very slowly and essentially in phase with each other. The 
contribution from chain 2 is exactly the same, and we find 



S» = 



1 - — \ Jd 2 x [(d x m) 2 +A1 2 



Inter-chain interaction is similarly evaluated 

2 

m , 



S± = 2Ja s 2 1 



d 2 xl 2 . 



(4.19) 



(4.20) 



Collecting terms thus obtained and integrating the 
field I out, we have S = J d 2 xC, where 



C 



1 

2fJ 



v(d\m) 



-(d 2 m) 



-TT^iivm ■ (dnm x d v m), 

87T 



(4.21) 
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where 



V. SUMMARY 



= 2iris = tti 
m 

s(m — 1) 
rn 



2Jaos 



to — 1 



V1 + VJ2, 
2 v/iT772, 



(4.22) 



with r = 1/J (= J\\/J±). 

It should be noted here that the topological term in 
the above action appears with the coefficient iri, which 
suggests that the present system should be massless, as 
we expected. Therefore, combining this result with the 
analyses in sees. II and III, we end up with the con- 
clusion that the present depleted system with spin sin- 
glet ground state should have massless spin excitations, 
belonging to the same universality class as the uniform 
s = 1/2 Heisenberg chain. 

Before concluding this section, a brief comment on the 
case of larger to is in order. As we mentioned, our present 
approach is valid for the case of rather small to. If we 
nai vely set to — > oo, bulk quantities such as g and v in 
eqr( 4.22 ) reduce to the formula derived by DeH'Aringa et. 



al.cj for the uniform non-depleted spin ladders, while the 
coefficient 8 remains 8 = iri. This result seems some- 
what paradoxical at first sight. However, one should 
note that the continuum limit in our approach means 
M = N/m — > oo. Therefore, however large to may 
be, the system always has the periodic array of impu- 
rities with finite density (infinite number of impurities 
in the thermodynamic limit), which may form the mass- 
less mode inside of the Haldane gap even for the limit of 
large to. This massless mode should be smoothly extrap- 
olated from the massless mode obtained for the small to 
case. In this sense, the limit to — > oo contains a subtle 
problem. We think that in the present approximation 
the topological term 8 — iri correctly reflects this mass- 
less mode whereas the values of g and v are determined 
by the host part with the Haldane gap in the large to 
limit. Therefore, to apply our analysis to the low-energy 
physics in the large to case consistently, we need to im- 
prove systematically the evaluation of the bulk quantities 
of g and v. In fact, we have checked in a preliminary 
calculation that the estimated values of the above bulk 
quantities are improved if we include not only the lowest 
spin mode but also higher modes. From these discus- 
sions, we believe that the present conclusion drawn from 
the above treatment applies also to the case of larger to at 
least qualitatively, although the approximation becomes 
worse when to becomes large. To improve our results in 
the large to case is an open problem to be explored in the 
future study. Note that in the case of larger to, an ap- 
propriate scheme to treat low-energy excitations has been 
considered recently by Nagaosa, et al. in the context of a 
randomly depleted spin laddercl. 



We have considered a class of the models for regularly 
depleted two-leg spin ladder systems, and investigated 
their low-energy properties. It has been proved rigor- 
ously that by the depletion this special class of the mod- 
els generates a new state which is characterized either by 
the presence of gapless excitations or a degenerate ground 
state. To see whether the massless states are really pro- 
duced by the depletion, we have investigated low-energy 
properties in terms of two different field theoretical ap- 
proaches. 

We have first proposed a scheme to describe low-energy 
properties by applying renormalization group methods 
to an effective weak-coupling model of Hubbard type, 
in which the effect of the depletion is incorporated in 
the strong coupling limit of the on-site interaction. This 
model has been analyzed in a weak coupling approxima- 
tion based on the one-loop RG method. Although we 
have seen how the tendency to a massless state is devel- 
oped when the effect of the depletion becomes strong, it 
has been found that our model is still massive for any 
finite couplings. This implies that it is not easy to de- 
scribe the effect of depletion in a naive continuum limit 
by a conventional method. 

To clarify this point, we studied the depleted spin 
model by mapping it to the non-linear sigma model. We 
have found via the spin-wave analysis that the dispersion 
relation of the lowest spin wave mode has indeed a linear 
spectrum in the case of the singlet ground state, which al- 
lowed us to mapping to the non-linear sigma model with 
a topological term. It has been shown that the coeffi- 
cient of this topological term is iri, which coincides with 
that of the models which have massless spin excitations 
such as the spin 1/2 Heisenberg chain. By combining all 
the analyses in sees. II, III and IV, we thus conclude that 
the periodic depletion of the present spin ladder systems 
should produce massless spin excitations. 

Finally we note that the drastic change from the spin- 
gap state to the gapless state observed in the present 
model reflects that the phase coherence in the wave func- 
tion is quite sensitive to the depletion. In this sense, 
the formation of the gapless spin state in our periodic 
model have the same basic origin as that with random 
non-magnetic impurities, although the low-temperature 
thermodynamics of the two is rather different u 
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APPENDIX A: CONTINUUM LIMIT OF 
INTERACTIONS 

In this Ap pen dix, w e d erive the continuum limit of 
interactions and (3.8). First, the continuum limit 

of the intra-chain interaction is calculated as 



H, 



(i) 



dx(ua + U'aOf) (2J? L J? R - 2Jf L J? R 



+e ilk - x 4 L ^ lR] 4 L ^ lRi + h.c.) (Al) 



Note here that 9j defined by eg. (2.2) can be rewritten as 



1 m 
m L — ' 



1 = 1 



,2-nilj jra 



(A2) 



Therefore, various kinds of terms behave as 

(ossillating terms), 



e 4ik F x _ e 2iri(l+l/m)j 



1 

m 



6j — — + (ossillating terms), 



)' x — + (ossillating terms), (A3) 



where i in 9^ is i = 1,2. Consequently, if we neglect 
oscillating terms, we have eq.([T|). 

Next consider the inter-chain coupling. By using the 
continuum limit of the spin operator, 



Sff'/oo -> J? L + J? R + (e 2 ^*N? 
where 



h.c. 



the inter-chain coupling terms are calculated as 



(A4) 



(A5) 



dx 1 



3 (i) 



a(2) 
3 



f (2k f oscillating terms) 
f (4kp oscillating terms) 



h.c. 



Note that 



1 - 



1 - 



3 (2) 



m/2 



711 * J 



t 4:7rilj jra 



(A6) 



(A7) 



Therefore, we have formulae similar to (A3) 

2 
m 



M")M 2 ') 
(i-^>)(i-<f) 

M")M 2) ) 



l 



(oscillating terms) , 



1 - 6»j 2) J e 2ifcF2; = (oscillating terms), 



(oscillating terms). 



(A8) 



By expressing N^N^ by L, 

N?N? + h.c. = \(L\L^ + h.c.) - \(L\l)« + h.c), 

(A9) 



we end up with eq.(3.8). 



APPENDIX B: OPERATOR PRODUCT 
EXPANSION 

The basic operator product expansion (OPE) for Fermi 
fields is 



SijSa/3 

2ir(z — w) ' 



(Bl) 



where z — vpT + ix. Here and in what follows, we neglect 
regular terms. Operator products should be normal- 
ordered, though we do not explicitly indicate this. Simi- 
lar formulae hold for right-moving currents, by replacing 
z -► z. Define f^ x and d^ x as [cr^/2, a v /2] = f^ x a x /2 
and o»o v = d^ x cr x . Namely, f ^ = = f^ = 



0, / 



a be 



,abc 



and d^ u = cP " 



d ab § c0 + it abc_ Then we haye 



8ir 2 (z-w) 2 



-L ^ 

27r(2 — w) 



8tt 2 (z — w) 2 

d^ x J x L {w) - d^ x J x L (w) 



M lL {z)M jL (w) 

J? L (z)Ll(w) ~ 
J* (*)2#M ~ 



+- 
(1 



27r(z — w) 



-1 



- 7° 



A.ir 2 {z — w) 2 ir(z — w) 
{S tl d^ x - 5 l2 d v ^)L x {w) 



J? L (z)M jL (w) ~ % 



47r(2: — u>) 
(1-2), 



2tt(z — w) 



(B2) 



The OPEs among various 4- Fermi interactions in the text 
follow from the above formulae. Using these OPEs, we 
can derive the one- loop order RG equation as follows: 
Suppose 



(B3) 



where the operators Oi represent the various 4-Fermi in- 
teractions with dimension 2 given above. The RG equa- 
tions for such perturbations are given by 

= (2 - x k )g k + 2vr 2 Cijkgtgj, (B4) 

i,3 
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where Xk is the dimension of the operator Ok and Cijk 
is the OPE coefficients defined by 



Oi(z)Oj{w) 



a 



^— ' \z — w\ 2 

k 1 



O k (w) 



APPENDIX C: BOSONIZATION 

The Fermi fields are bosonized as 
1 



IpiLa 
IpiRa 



'2-KCL 
1 

'27ra 



(B5) 



(CI) 



Coupling constants 



9-u 
9* 

Qxa 

{gt* - gt P )/2 
-{gta + gtp) 
Zgta 



Corresponding operators 



■} $ sin s/Airtpc sin y/4ivip c 

COS V^TYipf COS V^TVipJ 

cos ^/47r( ( o+ cos V^nOs 
cos y/4^6~ cos v^V^T 
cos VInpJ cos v^V 5 c 
COS ^4^V9+ cos \/47r¥>c 



Define 



TABLE II. Bosonized operators 



Via — V'iLa + Vi_Ra , 



and 

<Pi. 



1 

71 

1 
71 



(<ph - <pu), 



1 
1 

71 



(0»t + ^u)i 



(6» iT - iX ) 



(C2) 



(C3) 



For ladder systems, it is convenient to define furthermore 
the in- and out-of-phase Bose fields, 



(C4) 



Then, except for the gradient terms, we have the 
bosonized form of the various interactions in the text, 
summarized in the Table ||. 



APPENDIX D: SPIN WAVE HAMILTONIAN 

Here we su mm arize the expressions for the matrices h 
and A in eq.(4.4). Let us restrict ourselves to the J = 1 
case. First, the matrix h is diagonal and given by 



hij = -Sij x (# of the nearest neighbour for ith site), 

(Dl) 

wh ere zth site means its simplified notation used in 
eq.(4.4). For example, m — 6 case, we have 



h = ~diag(0, 2, 3, 2, 3, 2, 2, 3, 2, 0, 2, 3), (D2) 



and m = 8 case 



ft = ^diagtO, 2, 3, 3, 2, 3, 3, 2, 2, 3, 3, 2, 0, 2, 3, 3). (D3) 

Interaction part A is constructed as follows: First de- 
note it by m x m submatrices 



A= l S 1 

4 V lm <5 



(D4) 



where the matrix (5 represents intra-chain coupling, de- 
fined by 



/0 7 

7 7 
7 7 
7 



\7 



7\ 



7 
7 0/ 



(D5) 



with 



7 = exp(2irik/N), 



(D6) 



and l m is the m x m unit matrix, corresponding to the 
inter-chain coupling. The matrix A thus defined cor- 
responds to that of spin ladders without depletion. To 
include the effects of the depletion, the non-zero elements 
in the first row and column, as well as in the 3m/2 + 1th 
row and column should be set to 0. 
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